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Abstract 

We calculate the inclusive spectra of secondaries produced in soft (minimum bias) 
p+Pb collisions in the framework of Quark-Gluon String Model at LHC energy, and by 
taking into account the inelastic screening corrections (percolation effects). The role of 
these effects is expected to be very large at very high energies, and they should decrease 
the spectra about 3 times in the midrapidity region and increase them about 2 times in 
the fragmentation region at the energy of LHC. 

PACS. 25.75.Dw Particle and resonance production 
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1 Introduction 



The detailed investigation of p+Pb interactions at the highest LHC energy of \/s = 5 
TeV (4 TeV proton beam and 1.57 ■ A TeV Pb beam) makes part of the nearest plans 
by LHC. The investigation of soft p+Pb interactions is very interesting because it can 
give the final answer to the problem of inelastic shadow corrections [U |2] for inclusive 
particle production. 

In the present paper we present the predictions of the Quark-Gluon String Model 
(QGSM) [2111] for the rapidity and xp distributions of secondaries produced in p+Pb at 
y/s = 5 TeV. The QGSM quantitativey describes many features of high energy production 
processes, including the inclusive spectra of different secondary hadrons produced in 
high energy hadron-nucleon [SI El El IB] , hadron- nucleus [HJ [TD] , and nucleus- nucleus [UJ 
collisions. The model parameters were fixed by comparison of the theoretical calculations 
to the experimental data. 

In the case of interaction with a nuclear target the Multiple Scattering Theory 
(Gribov-Glauber Theory) is used, what allows to consider the interaction with the nu- 
clear target as the superposition of interactions with different numbers of target nucleons. 
However, it was shown in [TJ [2] that the description of the inclusive spectra of secon- 
daries produced in d+Au collisions at yfs = 200 GeV (RHIC) requires to account for 
the inelastic shadow corrections. These corrections are connected to the multipomeron 
interactions and they lead to the saturation of the inclusive density of secondary hadrons 
in the soft (low Pt) region, where the methods based on perturbative QCD cannot be 
used. The effects of the inelastic shadow corrections should increase with the initial en- 
ergy. The difference in the calculations with and without these effects at LHC energies 
is of about 3 times in the midrapidity region and of about 2 times in the fragmentation 
region. 

Other predictions for secondary production in p+Pb collisions at LHC energies (mainly 
for hard collisions can be found in [12j [131 HH US] • 

2 Inclusive spectra of secondary hadrons in the 
Quark-Gluon String Model 

In order to produce quantitative predictions a model for multiparticle production is 
needed. For that purpose we have used the QGSM [31 H] in the numerical calculations 
presented below. 

In the QGSM high energy hadron-nucleon and hadron-nucleus interactions are con- 
sidered as proceeding via the exchange of one or several Pomerons, and all elastic and 
inelastic processes result from cutting through or between Pomerons [16J. Each Pomeron 
corresponds to a cylinder diagram (see Fig. la) that, when cutted, produces two show- 
ers of secondaries as it is shown in Fig. lb. The inclusive spectrum of secondaries is 
determined by the convolution of diquark, valence quark, and sea quark distributions, 
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u(x,n), in the incident particles, with the fragmentation functions, G(z), of quarks and 
diquarks into secondary hadrons. Both functions u(x,n) and G(z) are determined by 
the corresponding Reggeon diagrams [IT] . The diquark and quark distribution functions 
depend on the number n of cut Pomerons in the considered diagram. 




Figure 1: a) Cylinder diagram (cylinder is shown by dash-dotted curves) corresponding to the one- 
Pomeron exchange contribution to elastic pp scattering, b) and the corresponding cut diagram which 
represents its contribution to the inelastic pp cross section. Quarks are shown by solid curves and string 
junctions by dashed lines. 



In the case of a nucleon target the inclusive spectrum of a secondary hadron h has 
the form [3]: 



oo 



dn 1 da x E da ^ , 

dy (Tinei dy a ind dx F ^ 

where the functions determine the contribution of diagrams with n cut Pomerons 

and w n is the probability for this process to occur [18]. Here we neglect the diffraction 
dissociation contributions which are important mainly for the secondary production in 
the large xp region that is not significant in the present calculations. 
For pp collisions 

<£{x) = fi q (x + ,n)-f^x^n) + f^x +} n)-f^ q (x. } n) + 2(n-l)f^x + ,n)-f^x.,n) , (2) 



x± = - Am^/ s + x 2 ± x] , (3) 

where f qq , f q , and f s correspond to the contributions of diquarks, valence quarks, and 
sea quarks, respectively. 

These contributions are determined by the convolution of the diquark and quark 
distributions with the fragmentation functions, e.g., 

f£{x + ,n)= [ u q (x u n) ■G ! q i (x + /x 1 )dx 1 . (4) 



3 



Figure 2: The distribution of probabilities to cut n Pomerons, w n in pp interactions (left), and the 
distribution of probabilities for the inelastic interaction of a proton with v lead nucleons, W p pb(v), in 
p+Pb interactions (right), at three initial energies, y/s = 5 TeV (solid curves), 200 GeV (dashed curves), 
and 20 GeV (dash-dotted curves). 



The diquark and quark distributions, as well as the fragmentation functions, are deter- 
mined by Regge intercepts [17]. The numerical values of the model parameters were 
presented in ref. [6]. 

The probabilities w n in Eq. (1) are the ratios of the cross sections corresponding to 
n cut Pomerons, a^ n \ to the total non-diffractive inelastic pp cross section, a n d [18J. 

The contribution of multipomeron exchanges in high energy pp interactions results in 
a broad distribution of w n . These distributions at three different energies y/s = 5 TeV, 
200 GeV, and 20 GeV are presented in the left panel of Fig. 2. 

In the case of nuclear target (or colliding beam) a projectile proton can interact 
with several nucleons. Let W p a{v) be the probability for the inelastic interactions of 
the proton with v nucleons and <Jpf od the total cross section of secondary production in 
a p+A collision. The calculated distributions of W p a{v) for proton- lead interactions at 
three different energies are presented in the right panel of Fig. 2. They were calculated 
in the framework of the Multiple Scattering Theory as 




(5) 



where 





coincides [HI |20j [2TJ [22] with the optical model expression [23] , and 




(7) 



with T(b) being the profile function of the nuclear target: 

/oo 
dz ■ p(b, z) , (8) 
-oo 

where p(r = y/b 2 + z 2 ) is the one-particle nuclear density. 
The average value of v has the well-known form: 

A ■ n pp 

<"> = ■ (9) 

prod 

We use the values a PP el = 72 mb and o v ^ d = 1900 mb at y/s = 5 TeV, so the numerical 
value of {v) p Pb turns out to be of about 7.9. 

In the calculation of the inclusive spectra of secondaries produced in pA collisions we 
should consider the possibility of one or several Pomeron cuts in each of the v blobs of 
proton-nucleon inelastic interactions. For example, in Fig. 3 it is shown one of the dia- 
grams contributing to the inelastic interaction of a beam proton with two lead nucleons. 
In the blob of the proton- nucleonl interaction one Pomeron is cut, and in the blob of the 
proton- nucleon2 interaction two Pomerons are cut. It is essential to take into account all 
digrams with every possible Pomeron configuration and its permutations. The diquark 
and quark distributions and the fragmentation functions here are the same as in the case 
of the interaction with one nucleon. 
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Figure 3: One of the diagrams contributing to the inelastic interaction of an incident proton with two 
target nucleons N\ and N2 in a pA collision. 



In particular, the contribution of the diagram in Fig. 3 to the inclusive spectrum is 
^-'T~ = 2-«-< 1 -< 2 -{/i(x +1 3)./:(x-,l)+ (10) 

a prod aXF 

+ f h q {x + , 3) • /£(*_, 1) + 3) • [/* 2) + 2) + 

+ 2./*(x_,2)]} . 
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The process shown in Fig. 3 satisfies [TS1 [2H [2E 122] the condition that the absorp- 
tive parts of the hadron-nucleus amplitude are determined by the combination of the 
absorptive parts of the hadron-nucleon amplitudes. 

Sometimes, in the case of hadron-nucleus collisions the values of xf have to be rescaled 
[2U125] to account for the part of the initial energy that is used for nucleus desintegration, 
though this correction becomes negligibly small at LHC energies. 

3 Inclusive spectra in p+A collisions at very high 
energy and inelastic screening (percolation) effects 

The QGSM gives a reasonable description [21 El [2S] of the inclusive spectra of different 
secondaries produced in hadron-nucleus and deuteron-nucleus collisions at energies ^/snn 
= 14—30 GeV. Also the experimental data of secondaries produced in Pb+Pb collisions 
at v /s/vtv = 17 GeV [23 I2H] were reasonably described (except for some problems with 
charged kaons) by a superposition model in [TT] . 

At RHIC energies the situation drastically changed. The spectra of secondaries pro- 
duced in pp collisions can be rather well described, but the RHIC experimental data for 
Au+Au collisions [211 [30] give clear evidence of the inclusive density saturation effects 
which reduce the inclusive density about two times in the central (midrapidity) region 
when compared to the predictions based on the superposition picture [3H [321 S3] . This 
reduction can be explained by the inelastic screening corrections connected to multi- 
pomeron interactions [I]. The effect is very small for integrated cross sections (many of 
them are determined only by geometry), but it is very important [1] for the calculations 
of secondary multiplicities and inclusive densities at high energies. 

Following the estimations presented in reference PQ, the RHIC energies are just of 
the order of magnitude needed to observe this effect. The inelastic screening can make 
PQ the inclusive density to decrease in the midrapidity region about two times at RHIC 
energies and about three times at LHC energies, with respect to the calculations without 
inelastic screening. 

However, all estimations are model dependent. The numerical weight of the contribu- 
tion of the multipomeron diagrams is rather unclear due to the many unknown vertices 
in these diagrams. The number of unknown parameters can be reduced in some models, 
and for example in reference [1] the Schwimmer model [31] was used for the numerical 
estimations. 

Another (model dependent) possibility to estimate the contribution of the diagrams 
with Pomeron interaction comes [351 ESI EH SHI ESJ from percolation theory. The perco- 
lation approach and its previous version, the String Fusion Model [101 El H2] predicted 
the multiplicity suppression seen at RHIC energies, long before any RHIC data were 
taken. 

The percolation model also provides a reasonable description of the transverse mo- 
mentum distribution (at low and intermediate px) including the Cronin effect and the 
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behaviour of the baryon/meson ratio [HH EH US]- Most of the effects predicted by per- 
colation can be seen as a direct consequence of the strong colour fields produced in the 
collision. This feature is common to other approaches as the Colour Glass Conden- 
sate [46j H7], where a pr scaling is also obtained [48J. 

In the percolation approach one assumes that if two or several Pomerons overlap in 
transverse space, they fuse in only one Pomeron. When all quark-gluon strings (cut 
Pomerons) are overlapping, the inclusive density saturates, reaching its maximal value 
at a given impact parameter. This approach has only one free parameter rj, called the 
percolation parameter 

V = N s .^.{r(y)), (11) 

with N s the number of produced strings, r s the string transverse radius, and R the radius 
of the overlapping area. The factor (r(y)) accounts for the fact that the parton density 
near the ends of the string is smaller that in the central region, where we fix r(0) = 1. 
At large rapidities we have N s strings with differen parton densities, rj(y), and 

N s 

N s -(r(y)) = J2n(y) • (12) 

i=l 

As a result, the bare inclusive density dn/dy\ bare gets reduced, and we obtain 

dn/dy = F{rj) ■ dn/dy\ bare , (13) 

with [49] 

F(V) = \J^~ ■ ( 14 ) 

In order to account for the percolation effects, it is technically more simple [2] to 
consider in the central region the maximal number of Pomerons that can be emitted 
by one nucleon, n max . By doing this, all model calculations become rather simple since 
above the critical value n max every additional Pomeron cannot contribute to the inclusive 
spectrum (it fuses with some another Pomeron). In this scenario we obtain a reasonable 
agreement with the experimental data on the inclusive spectra of secondaries at RHIC 
energy (see [2] with n max = 13). In the following calculations we use the simplest version 
of this approach, in which we neglect the energy dependence of n max , i.e. we use a fixed 
n ma x = 13 at LHC energy (see discussion in the section Conclusions). 

This scheme seems closer to the point of view of the Parton Model [50l I5T] . the 
numerical difference between the present secheme and [501 [51] coming from the fact that 
the ratio r 2 / R 2 in Eq. (11) is rather small [2], so the percolation parameter, 7], is also small 
and many n < n max independent Pomerons can exist before percolation plugs in. That 
explains why the effects of Pomeron (secondary particle) density saturation are small at 
fixed target energies and they become visible only starting from RHIC energies. The 
similar energy dependence in the Reggeon calculus pQ is connected to the suppression 
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of the diagrams with multipomeron interactions due to the nuclear longitudinal form 
factor. 

In the following calculations one additional effect is also accounted for, namely the 
transfer of the baryon charge to large distances in rapidity distances through the string 
junction effect [7J [52j |53j [Ml ESI [561 EZ]- This effect leads to an asymmetry in the 
production of baryons and antibaryons in the central region, that is non-zero even at 
LHC energies (SB]- In the calculations of these effects we use the following values of the 
parameters [56J: 

a SJ = 0.5 and e = 0.0757. (15) 



4 Rapidity spectra of secondaries at LHC energies 

In the lead-nucleus fragmentation region the contributions of intranuclear cascade and/or 
Fermi- motion of bound nucleons exist. To avoid these contributions we consider only 
the central (midrapidity) and the proton beam fragmentation regions y cm (NN) > 0. 

The rapidity distributions of all charged secondaries dn c h/dy cm , where y cm is deter- 
mined in the frame of a incident proton interaction with one lead nucleon, calculated in 
the QGSM are presented on the left panel of Fig. 4. The solid curve shows the results 
obtained by accounting for the inelastic screening (percolation) corrections, the dashed 
curve shows the results obtained without these corrections, and the QGSM predictions 
for pp collisions at ^/snn = 5 TeV are presented by the dotted curve. 

In the right panel of Fig. 4 we present the predictions for the nuclear modification 
factor 

R h (y ) ~ dUch/ ^ cm ( pPb ) (16) 
Ck Cm dn ch /dy cm (pp) 

Here again, the solid curve shows the QGSM predictions when accounting for the inelastic 
screening and the dashed curve corresponds to the results without these corrections. 

In Fig. 4 one can clearly see that the values of dn c h/dy, as well as those of the ratios 
Rch{pPb/pp), differ in the midrapidity region for about 3 times between the calcula- 
tions with and without percolation efects, what is in agreement with the predictions 
in reference [TJ. These differences decrease with rapidity, and at large y cm , i.e. in the 
fragmentation region, these differences change sign, the values R c h{pPb/pp) becoming 
smaller than unity in agreement with the nuclear shadowing in the inclusive spectra 

[aiEEniEH]. 

The effect of the nuclear shadowing at large rapidities differs for about two times 
for the calculations with and without percolation efects, what is a direct consequence of 
energy conservation. The behaviour of the spectra in the proton fragmentation region is 
discussed in more detail in the next section. 

The results of similar calculations for secondary protons are shown in Fig. 5. Here 
again, the difference in the predictions in the midrapidity region for the calculations with 
and without percolation efects is of about 3 times, and these differences decrease with 
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Figure 4: Rapidity distributions (left panel) and nuclear modification factors (right panel) for all 
charged secondaries produced in p+Pb collisions at ^/smn — 5 TeV. Dashed curves show the QGSM 
predictions without inelastic screening (percolation) corrections, solid curves show the QGSM predictions 
with these corrections, and the dotted curve in the left panel shows the QGSM predictions for pp collisions 
at y/s NN — 5 TeV. 

rapidity. The rather complicate behaviour of the spectra in the fragmentation region is 
also discussed in the next section. 

The calculated results of dn/dy cm distributions for secondary A and A are shown in 
the upper panels of Fig. 6. The nuclear modification factors R\(pPb/pp) and Rx{pPb/pp) 
practically coincide inside the 10—15% accuracy margin, and they are presented in the 
lower panel of Fig. 6. 

As a result, we can predict very similar behaviour of the nuclear modification factors 
for different secondaries in the midrapidity region and similar dependences of these fac- 
tors on y cm . At the same time, these factors are about 3 times different when they are 
calculated with and without inelastic screening corrections. 

5 Feynman-x spectra at LHC energies 

The xf variable is mainly suited for the consideration of the inclusive spectra in the frag- 
nentation region. For the case of the distributions [x E joi ne i) • {da/dxp), xf is determined 
in the cm. frame of the interaction of an incident proton with one lead nucleon. 

The distributions (x E /a ine i) ■ (da/dxp) calculated in the QGSM are presented in 
the left panel of Fig. 7. The solid curve shows the result obtained with accounting for 
the inelastic screening (percolation) corrections, the dashed curve represents the result 
without these corrections, and the QGSM predictions for pp collisions at v /% 5 " = 5 TeV 
are presented by the dotted curve. 
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Figure 5: Rapidity distributions (left panel) and nuclear modification factors (right panel) for sec- 
ondary protons produced in p+Pb collisions at ^/sjvjv = 5 TeV. The correspondence of the curves is the 
same as in Fig. 4. 



In the right panel of Fig. 7 we present the predictions for the nuclear modification 
factor 

, x _ (l/Oprod) • (d<T ch /dXp) p pb 
(1/cTmeO ■ (d(Tch/dx F ) pp 

Here again, the solid curve shows the QGSM predictions with accounting for the inelastic 
screening and dashed curve shows the result without these corrections. 

The results of the corresponding calculations for the secondary protons are shown 
in Fig. 8, where the effects of nuclear shadowing in the fragmentation region are very 
clearly seen, being about two times different at xf > 1/2 for the calculations with and 
without percolation efects. 

Some predictions of the nuclear shadowing in the fragmentation region were ob- 
tained [5H] in the framework of the Additive Quark Model (AQM). In this model the fast 
baryon is considered as formed by three constituent quarks, each carrying out around 
1/3 of the total baryon momentum (sea quarks and gluons are assumed to be effectively 
included inside the constituent quarks). Let a secondary baryon with xf ~ 2/3 to be 
formed by two constituent quarks from the initial baryon, the examples being the re- 
actions p — > p, p — > A, etc. The simplest mechanism in the framework of AQM to get 
this is the recombination of two quark-spectators of the incident baryon with one newly 
produced quark from the sea. 

The A-dependence of such a processes is determined by the probability to have two 
quark-spectators in pA collisions. If the constituent quarks interact with the target 
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Figure 6: Rapidity distributions for secondary A (left up) and A (right up) produced in p+Pb col- 
lisions at y/sNN — 5 TeV, and the corresponding nuclear modification factors for them (low). The 
correspondence of the curves is the same as in Fig. 4. 



independently, the cross section of their inelastic interaction with a proton target is 



qN 
inel 



NN 
'•'inel i 



(18) 



and the probability that one constituent quark will interact with the target nucleus while 
two other constituent quarks are spectators can be written as: 



V VA = _L_.J d 2 b . e -2^-T( 6 ) . _ e -<- .T( 6 ) 
^ 'prod 



(19) 



where T(b) is determined by Eq. (8). 
The AQM predicts an A dependence 



R 



.(x F ~2/3) = Vr A , 



(20) 
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Figure 7: Feynman-x distributions (left panel) and nuclear modification factors R c h{xf) (right panel) 
for all charged secondaries produced in p+Pb collisions at ^smn = 5 TeV. The correspondence of the 
curves is the same as in Fig. 4. 

and such a simple relation is actually in good agreement [TDJ EO] with the experimental 
data. 

The AQM prediction from Eq. (20) is shown in the right panel of Fig. 8 by a dash- 
dotted straight line, and it is in agreement with the QGSM calculation without inelastic 
shadow corrections. Shadow corrections (elastic as well as inelastic ones) also exist in 
the AQM but they are not accounted for in eqs. (18) and (19). 

The calculated results for the (x e I '&inei) • {do a/ 'dx f) distributions, and the correspond- 
ing nuclear modification factors R\{xp), are shown in Fig. 9, where the AQM prediction 
from Eq. (20) is presented in the right panel of Fig. 9 by the dash-dotted straight line. 

In Fig. 9 one can clearly appreciate that the spectra of A at large Xp are very small. 



6 Conclusions 

The percolation (inelastic screening) effects are too small to be observed at fixed target 
energies ^/snn = 15—40 GeV by comparing the model calculations to experimental data. 
At RHIC energies these effects in dn c h/dy(y = 0) become to be a factor 1.5—1.7 [2], and 
they can increase until a factor 3 at LHC energy. A so large effect at the LHC energies is 
connected to the fact that at ^stvtv = 5 TeV, (tinn) ~ 3 and (vNPb) ~ 8, so the average 
value of Pomerons (~ 24) in the minimum bias p+Pb collision is significantly larger than 
n max = 13. The possible weak energy dependence [2] of n max = 13 cannot significantly 
change our results. 

To prove the existence of the inelastic screening effects it is enough to measure the 
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Figure 8: Feynman-x distributions (left panel) and nuclear modification factors (right panel) for 
secondary protons produced in p+Pb collisions at ^snn = 5 TeV. The correspondence of the curves is 
the same as in Fig. 4. The AQM predictions (eqs. (19) and (20)) are shown by the dash-dotted straight 
line. 



nuclear modification factors R(y C m) of Eq. (16) in the midrapidity region and the ratio 
R(xf) of Eq. (17) in the fragmentation region. The first one should be several (3 or 4) 
times smaller than the value (y) given by Eq. (5), and the second one atXi?~l/2 — 3/4 
should be about 2 times larger than the QGSM or AQM predictions. The experimental 
confirmation of the inelastic screening (percolation) effects for the inclusive spectra can 
be considered as a natural "bridge" between the mechanisms of the hadron density 
saturation in soft and hard processes. Our results are in reasonable numerical agreement 
with those in ref. [6T| . 

Other scenarios in the percolation or in the Colour Glass Condensate approaches (see 
[2] and references therein) exist, in which the saturation effects have a different energy 
dependence, so the experimental investigation of inclusive spectra in p+Pb collisions at 
LHC can unravel the question of which among the possible scenarios are excluded or 
confirmed by the experiment. 
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Figure 9: Feynman-a; distributions (left) and nuclear modification factors (right) for secondary A 
produced in p+Pb collisions at y/s~NN = 5 TeV. The correspondence of the curves is the same as in 
Fig. 4. The AQM predictions from eqs. (19) and (20) are shown by a dash-dotted straight line. 
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